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Net-Charge fluctuations in a hadron gas are studied using an effective liadronic interaction. Tlie 
empfiasis of tliis work is to investigate tfie corrections of fiadronic interactions to the charge fluctua- 
tions of a non-interacting resonance gas. Several methods, such as loop, density and virial expansions 
are employed. The calculations are also extended to SU{3) and some resummation schemes are con- 
- - - , sidered. Although the various corrections are sizable individually, they cancel to a large extent. As 

. a consequence we find that charge fluctuations are rather well described by the free resonance gas. 
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The study of cvcnt-by-cvcnt fluctuations or more generally fluctuations and correlations in heavy ion collisions has 
recently received considerable interest. Fluctuations of multiplicities and their ratios [1], transverse momentum [2-5] 
and net charge fluctuations [6-9] have been measured. Also flrst direct measurements of two particle correlations have 
been carried out [10, 11]. 

, Conceptually, fluctuations may reveal evidence of possible phase transitions and, more generally, provide information 
^\ • about the response functions of the system [12]. For example, it is expected that near the QCD critical point long 
^ range correlation will reveal themselves in enhanced fluctuations of the transverse momentum (pt) per particle [13]. 
■ Also, it has been shown that the fluctuations of the net charge are sensitive to the fractional charges of the quarks in 
O ; the Quark Gluon Plasma (QGP) [14, 15]. 

(~| . Most fluctuation measures investigated so far are integrated ones, in the sense that they are related to integrals 
of many particle distributions [16]. Examples are: Multiplicity, charge and momentum fluctuations which are all 
related to two-particle distributions. These integrated measures have the advantage that they can be related to well 
defined quantities in a thermal system. For example, fluctuations of the net charge are directly related to the charge 
^ ] susceptibility. However, in an actual experiment additional, dynamical, i.e non-thermal correlations may be present 
which make a direct comparison with theory rather difficult. This is particularly the case for fluctuations of the 
transverse momentum, where the appearance of jet like structures provides nontrivial correlations [10, 17, 18]. These 
need to be understood and eliminated from the analysis before fluctuation measurements can reveal insight into the 
$-H ' matter itself. 

In this article we will not be concerned with the comparison with experimental data, and the difficulties associated 
with it. We rather want to investigate to which extent interactions affect fluctuations. Specifically, we will study the 
fluctuations of the net electric charge of the system, the so-called charge fluctuations (CF). CF have been proposed 
as a signature for the formation of the Quark Gluon Plasma (QGP) in heavy ion collisions [14, 15]. Refs. [14, 15] 
note that CF per degree of freedom should be smaller in a QGP as compared to a hadron gas because the fractional 
charges of the quarks enter in square in the CF. Using noninteracting hadrons and quarks, gluons, respectively, it 
was found that the CF per entropy are about a factor of 3 larger in a hadron gas than in a QGP. The net CF per 
entropy has in the meantime been measured [6-9] . At RHIC energies the data are consistent with the expectations of 
a hadron gas, but certainly not with that of a QGP. This might be due to limited acceptance as discussed in [19-21]. 

The original estimates of the net charge fluctuations per entropy in the hadron gas [14, 15] have been based on a 
system of noninteracting particles and resonances. While this model has been proven very successful in describing the 
measured single particle yields [22, 23], it is not obvious to which extent residual interactions among the hadronic states 
affect fluctuation observables. For example, in the QGP phase, lattice QCD calculations for the charge susceptibility 
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and entropy-density differ from the result for a simple weakly interacting QGP. Their ratio, however, agrees rather 
well with that of a noninteracting classical gas of quarks and gluons [12, 14, 24-26]. As far as the hadronic phase is 
concerned, lattice results for charge fluctuations are only available for systems with rather large pion masses [24, 26]. 
In this case, an appropriately rescaled hadron gas model seems to describe the lattice results reasonably well [27]. 
Lattice calculations with realistic pion masses, however, are not yet available. Thus, one has to rely on hadronic 
model calculations in order to assess the validity of the noninteracting hadron gas model for the description of CF. 
In Ref. [28] the electric screening mass m^j which is closely related to CF has been calculated up to ncxt-to-lcading 
(NLO) order in tttt interaction. However, the fact that thermal loops pick up energies in the resonance region of the 
TTTT amplitude where chiral perturbation theory is no longer valid leads to large theoretical uncertainties. 

It is the purpose of this paper to provide a rough estimate of the effect of interactions in the hadronic phase, in 
particular the effect of the coupling of the p-meson to the pions. Since p-mesons are strong resonances which carry 
the same quantum numbers as the CF this should provide a good estimate for the size of corrections to be expected 
from a complete calculation; the latter will most likely come from lattice QCD, once numerically feasible. 

As a first step we will consider the case of a heavy p-meson or, correspondingly, a low temperature approximation. 
In this case the p-meson is not dynamical and will not be part of the statistical ensemble. It will only induce an 
interaction among the pions which closely corresponds to the interaction from the lowest order (LO) chiral Lagrangian. 

Although the temperatures in the hadronic phase are well below the p-mass, it is interesting to estimate the residual 
TTTT correlations introduced when this resonance is treated dynamically. Special attention is paid to charge conservation 
and unitarity. In addition, we will investigate the importance of quantum statistics. Finally an extension to strange 
degrees of freedom is provided. 

This paper is organized as follows. After a brief review of the charge fluctuations we introduce our model Lagrangian 
and discuss the heavy rho limit. Next we discuss the treatment of dynamical p- mesons up to two- loop order and 
compare with the results obtained in the heavy rho limit. Then, the effect of quantum statistics and unitarity is 
discussed. Before we show our final results including strange degrees of freedom, we will briefiy comment on possible 
resummation schemes. 



II. CHARGE FLUCTUATIONS AND SUSCEPTIBILITIES 



Before turning to the model interaction employed in this work, let us first introduce some notation and recall the 
necessary formalism to calculate the CF (for details, see, e.g., Ref [12]). 

In this work we will consider a system in thermal equilibrium. In this case the charge fiuctuations {SQ'^') are given 
by the second derivative of the appropriate free energy F with respect to the charge chemical potential /i: 

m - = -VTxQ. (1) 

Here, T (V) is the temperature (volume) of the system and XQ is the charge susceptibility, which is often the preferred 
quantity to consider, particularly in the context of lattice QCD calculations. Equivalently, the CF or susceptibility 
are related to the electromagnetic current-current correlation function [29, 30] 

U^,{u;,k) = i J dtd^xe-'f"*-"^) (J,,(x,t)J,(0)) (2) 



((5g2) =FTnoo(u> = 0,k^O) ^FTm^i (3) 

which is illustrated for scalar QED in Appendix A. Relation (3) also establishes the connection between the CF and 
the electric screening mass md. 

As noted previously, the observable of interest is the ratio of CF over entropy 

* - W- C) 

Given a model Lagrangian, both CF and entropy can be evaluated using standard methods of thermal field theory 
(see e.g. [30]). CF are often evaluated via the current-current correlator using thermal Fcynman rules; evaluating the 
free energy and using relation (1) will lead to the same results as will be demonstrated in Sec. V. 
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Let us close this section by noting that in an actual experiment a direct measurement of the entropy is rather 
difficult. However, the number of charged particles (iVch) in the final state is a reasonable measure of the final state 
entropy. Therefore, the ratio 

has been proposed as a possible experimental observable for accessing the CF per degree of freedom. For details and 
corrections to be considered see Ref. [12] and references therein. In this article we will concentrate on the "theoretical" 
observable Ds defined in Eq. (4). 



III. MODEL LAGRANGIAN AND tttt INTERACTION IN THE HEAVY p LIMIT 



As already discussed in the Introduction, in this work we want to provide an estimate of the corrections to the CF 
introduced by interactions among the hadrons in the hadronic phase. Since it is impossible to account for all hadrons 
and their interactions, we will concentrate on a system of pions and p-mesons only, with some extensions to SU{2>) in 
later sections. A suitable effective Lagrangian for this investigation is the "hidden gauge" approach of Refs. [31, 32]. 
In this model the p-meson is introduced as a massive gauge field. The Trp interaction results from the covariant 

[p^, $] acting on the pion field U{x) = exp[i$(a;)//7r] in the LO chiral Lagrangian 



derivative Z?^$ 



2 



- ^ Tr [d^U^df'U + M(U + U^)] 



(6) 



by the replacement D^. Here, 



$ = 





(7) 



and = 93 MeV is the pion decay constant. An extension of the heavy gauge model to SU{?>) has been applied for 



vacuum and in-medium processes (see, e.g. 
in Sec. VII A. 

The resulting ■np interaction terms are 



Refs. [33, 34]) and is straightforward [35]. This extension is considered 



(8) 



and 



f^Tr([p„<I>]^ 



(9) 



Chiral corrections to the interaction in Eq. (8) are of 0{p^) or higher as pointed out in Ref. [36]. The interaction 
of Eq. (9) does not depend on the pion momentum, thus violating the low energy theorem of chiral symmetry [36]. 
Nevertheless, this term is required by the gauge invariance of the p- meson [37] and in fact cancels contributions in 
the pole term and crossed pole term of irp scattering via Eq. (8). 

To leading order in the pion field we, thus, have the following model Lagrangian: 



piTTT 



c 



ppTTTT 



(10) 



with the free field terms 



$2 



= -iTr(G^.Gn + ^Tr(m>^p^), 



(11) 



and the interaction terms Cpj^TT and Cppinr as given in Eq. (8) and (9), respectively. For the kinetic tensor of the p, 
Gpv = dp,pv — dvPp., we restrict ourselves to the Abelian part; a non-Abclian p would lead to additional Zp and 4p 
couplings. In the thermal loop expansion this would result in closed p loops which arc kincmatically suppressed. The 
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coupling constant g is fixed from the p — > tttt decay to be 17 = gp-n-n = 6 and and we use m-^ = 138 MeV and nip = 770 
MeV throughout this paper. 

As a first approximation, we start with the low temperature limit of the tttt interaction in which the p-mcson 
mediates the interaction of the pions but does not enter the heatbath as an explicit degree of freedom. To this 
end we construct an effective interaction based on s-, t-, and u-channel p-meson exchange as given by second order 
perturbation theory of the interaction Cp-mi. Furthermore, we assume that the momentum transfer of two pions 
interacting via a p is much smaller than the mass of the p-meson, >> fc^, i.e. we replace the propagator of the 
exchanged p-meson by -I /ml. Thus, we arrive at the following effective interaction 



Ct. = ir^{{^~ d, 2(7r" dp 7r+)(^° ^-)). (12) 



2 



2m'i 



p 

Note that in this limit, subsequently referred to as the "heavy p limit" the ppmr term from Eq. (9) does not contribute 
at order g". 

The effective Lagrangian of Eq. (12) shows the identical isospin and momentum structure as the kinetic term of 
Eq. (6) at 1//^. However, comparing the overall coefficient one arrives at 

ml = 3/Jg2 (13) 

which differs by a factor of 3/2 from the well known KSFR relation [38] = 2f^g^. We should point out the 
same factor has been observed in the context of the anomalous 7777777 interaction [39]. As discussed in more detail in 
Appendix B 1 the KSFR relation is recovered if one restricts the model to the s-channel diagrams for the isovector 
p-wave (Til) amplitude. Once also t- and w-channels are taken into account the factor 3/2 appears. For this study, 

we prefer the interaction (12) over £^r7r from Eq. (6) as it delivers a better data description at low energies in the 
p-channel (see Appendix B 1). The simplification from the "heavy p" limit of the 7777 interaction will later be relaxed 
in favor of dynamical p-exchange. However, the interaction in the heavy p limit will still serve as a benchmark for the 
more complex calculations. 

Since we are interested in the electromagnetic polarization tensor, the interaction of Eq. (12), together with the 
kinetic term of the pion, is gauged with the photon field by minimal substitution, leading to 

2 2 

+ (77-i:'„77+- 77+1); 77-)^-^ (77° i:i'^77+- 77+9'^ 77") (77° 15*77- -77-9^77°) (14) 

zmp mp 

with the covariant derivative of the photon field Dp = dp + ieAp, e > 0, and the photon field tensor F'^'^, leading to 
the 77777 and "/"fmr interactions of scalar QED, plus 777777777 and jjmrmT vertices. Vector meson dominance leads to 
7p'' mixing as pointed out, e.g., in Ref. [33], additionally to the vertices of Eq. (14). However, since the correlator 
of Eq. (3) is evaluated at the photon point, the form factor is unity and the process 7 ^ p° ^ 7777 which emerges 
in the systematic approach of Ref. [40] does not contribute to the 77777 coupling. Thus, no modification of Eq. 
(14) is required. Note also that the anomalous interaction providing 777P vertices [33] does not contribute in the 
long- wavelength limit studied here. This follows a general rule noted in Ref. [25]. 

In the following chapters, the p will be also treated dynamically. The interaction with the photon is then given 
by the scalar QED vertices from above, plus a ^pmr vertex which is obtained from Eq. (8) by minimal substitution. 
With the same procedure the direct 7p interaction is constructed from Eq. (11), leading to the vertices 





p-ia,k) p-{u,k') ^eik'^gp^ + k'^gp.-ik + k'rga.), P^W = 2e^ {gppga., - gp^g^^p) (15) 

in the imaginary time formalism. 



IV. CHARGE FLUCTUATIONS AT LOW TEMPERATURES 



Having introduced the effective interaction in the heavy p limit, we can evaluate the correction to the CF due to 
this interaction. Before discussing the results let us first remind the reader about the basic relations for CF in a 
noninteracting gas of pions and p-mesons. 
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FIG. 1: Photon selfenergy at for the free pion gas (C, D) and the free p gas {C, W). 



A. Charge fluctuations for free pions and p- mesons 

In order to illustrate the relations of Sec. II and to establish a baseline it is instructive to calculate Ds from Eq. 
(4) for the free pion gas in two ways: once via Eq. (3) and also directly from statistical mechanics. The interaction 
from Eq. (14) reduces to scalar QED in the zeroth order in g. To order the selfenergy is given by the set of gauge 
invariant diagrams in Fig. 1 and reads 

n"°(fc° = 0,k^O) =e2(C + L'), {5Q^) ^ e^TV {C + D) (16) 

according to Eq. (3) with 

oo oo 

C ^ \ f dp u; n[oj], D^^fdpp^^ (17) 



TT^ J TT 

b 



where uj = y/jp + rn^ the pion energy, n[uj] ~ l/(exp (/Jw) — 1) the Bose-Einstein factor, and /3 = l/T. The CF from 
Eq. (16) can also be derived from statistical mechanics, 

{SQ^)=e'T' -^\^=o\ogZ, (18) 



Both the photon selfenergy and Eq. (18) lead to the same CF also at the perturbative level as will be seen in Sec. V. 
The value for the chemical potential of fii = ±yU in Eq. (19) corresponds to charged pions and /Xi = is assigned to 
neutral pions which do not contribute to the CF but to the entropy S = d{T log Z)/dT of the free gas, 

oo 

So = ^^ fdpp' n[u;] ( 3lo + ^) . (20) 







In the high temperature limit, or for massless pions. the relevant thermodynamical quantities are given by 

(,0^> = £!It^ S = ?^T^ (A..,.) = H«2It3 (21, 

3 15 TT^ 

where (A^ch) is defined as in Ref. [14]. For the quantity Ds from Eq. (4) we obtain Ds ~ 0.185 for massive free pions 
at T = 170 MeV and Ds = 0.253 for massless pions. For Dc from Eq. (5), the values are 4.52 and 5.47, respectively. 

The classical (Boltzmann) limit is obtained by replacing the Bose-Einstein distribution n in Eqs. (17) and (20) by 
the Boltzmann distribution tt-b = exp(— /3w). In this case at T = 170 MeV we obtain Ds = 0.156 and Ds — 1/6 for 
massive and massless pions, respectively. For all masses and temperatures, Z?c = 4 the classical limit. For a QGP 
made out of massless quarks and gluons, Ds = 0.034, following the same arguments as in [14]. This is about a factor 
of five smaller than a pion gas. 

The CF for the free p gas are given by the diagrams with the double lines in Fig. 1. With the p propagator 

nt"" = i (gt^" - ^^^^ I (22) 
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(eff2) 



(eff3) 



(eff4) (effS) (eneff) 

FIG. 2: Selfenergy for tttt interaction in the heavy p limit at order e^g^ (diagrams (effl) to (effS)). Expansion of log^ at 
for the calculation of the entropy (diagram (eneff)). 

TABLE I: Static selfenergy n°°{k° = 0, k ^ 0) from Fig. 2 with C and D from Eq. (17). 



Diagram Contribution 



(eff2) -^^DiD-3C-P^^iC-D)) 

(eff3) +11^D{2D-C) 
mf, 

(eff4) CD 
mf, 



and the interaction from Eq. (15) the photon selfenergy turns out to be 

n°" (fc" = 0, k ^ 0) = [CP + DP] (23) 

where the upper index means that the pion mass is substituted by the p-mass in C and D from Eq. (17). The factor 
of three corresponds to the sum over the physical polarizations of the p. The same factor also appears in log Zo of 
Eq. (19) for the p. 



B. TTTT interaction in the heavy p limit to order e^g'^ 

At order e^g^ the Feynman rules derived from the heavy p limit Eq. (14) lead to the set of five diagrams (effl) to 
(eff5) depicted in Fig. 2. They are gauge invariant as shown in Appendix C4. The summation over Matsubara fre- 
quencies has been performed by a transformation into contour integrals following Ref. [30]. The limit (fc" = 0, k — > 0) 
for the external photon has to be taken before summation and integration, as discussed in Appendix A. The loop 
momenta factorize so that the diagrams of Fig. 2 can be expressed in terms of the quantities C and D from Eq. (17) 
as shown in Tab. L The sum of the diagrams is cast in a surprisingly simple form, 



^n,(fc° = 0,k^ 0) 



^{C + Df-PD^{C-D) 



(24) 



The entropy correction at is calculated from logZ given by diagram (eneff) in Fig. 2, 



^1 



3gW fm^ 
2T \ TO,, 



D{C + D). 



(25) 
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p + k 




(eff5) 

( ) — >• i it^ ] 

(enl) — ' (eneff) 

FIG. 3: Overview of the relevant two-loop diagrams at e^g^ for the photon selfenergy and at for the entropy. Diagram 
(la2) is calculated in detail in Appendix C, where also the results for all other diagrams and a proof of gauge invariance are 
found. The diagrams on the right hand side (effl-effS) correspond to the heavy p limit of the ones given on the left. This limit, 
indicated with arrows, is numerically shown in Appendix C 1. 

Note that using the LO chiral Lagrangian from Eq. (6) instead of the tttt interaction in the heavy p limit, results would 
simply change by a factor of (2/3)^, up to tiny corrections, which are due to higher order contributions involving the 
chiral symmetry breaking term ~ M from Eq. (6). Numerical results can be found in Sec. VI B, which supersede our 
findings from Rcf. [41]. 

V. THE p-MESON IN THE HEATBATH 

In this section we will relax the assumption of a heavy non-dynamical p-meson. This will allow for an estimate of 
the CF from the residual interactions of the p when this particle is treated as an explicit degree of freedom. It will 
also avoid some problems induced in the calculation from vertices of higher order in momenta such as encountered 
in the calculation in Ref. [28] (sec discussion in Sec. VIA, VI B). The pnn interaction from Eq. (8) involves 
vertices only linear in momentum and a smoother temperature dependence is expected. 

We start with the calculation of the diagrams in the first two columns of Fig. 3 because this subset corresponds 
to the heavy p limit from Sec. Ill; by increasing the p-mass from its physical value to infinity in these diagrams, the 
previous results from Tab. I are recovered as illustrated in Appendix C 1. Note that there is no need to include jp^ 
mixing or anomalous vertices as we have already seen in Sec. III. 
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X 
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-J 



(2c)'' 

(3c) 

(5c2^^...,-''' 

3E 

(7c) 



(8c) 



FIG. 4: Additional, sub-leading, diagrams at e^g^ with direct {■y)'ypp couplings and with pprnv interaction. Also, the diagrams 
which vanish are shown [(lc2), (2c), (3c), (4c), (5c2), (7c), (8c)]. 



Here and in the following sections, the p is treated as a stable particle (propagator from Eq. (22)) and we ignore 
imaginary parts at the cost of unitarity violations as will be discussed in Sec. VIA. A p with finite width would 
induce problems concerning gauge invariance: one would have to couple the photon to all intermediate p selfenergy 
diagrams that build up the p width in the Dyson-Schwinger summation. In principle, this is possible — see the last 
part of Sec. VIII — but goes beyond the scope of this work. The results for the diagrams with dynamical p from Fig. 
3 are found in Eq. (C2,C4) and Fig. 16 of Appendix C 1, together with a detailed calculation of one of the diagrams 
and a discussion of the infrared divergences. In Appendix C 4 the gauge invariance of the diagrams is shown. 

At order e^g^ there are additional diagrams with direct ^pp and 77pp couplings from Eq. (15) and also with the 
ppiTTT coupling from Eq. (9) which is required by the gauge invariance of the p-meson. The resulting diagrams are 
displayed in Fig. 4. Some of these diagrams contain more than one p-propagator. They are sub-dominant because 
every p propagator counts as l/m^. Furthermore, Fig. 4 shows diagrams which have a closed pion loop with only 
one vertex of the pmr type (see, e.g., diagram (2c)). The latter diagrams vanish due to the odd integrand in the loop 
integration. The set of diagrams from Figs. 3 and 4 is complete at order e^g^. 

The non-vanishing diagrams from Fig. 4 are best calculated by evaluating the corresponding partition function, 
\ogZ, at finite chemical potential p and differentiating with respect to p [25, 28] (see also Eq. (18)). For a calculation 
at finite p we first convince ourselves that for the simple interaction from Eq. (12) the use of Eq. (18) leads to the 
same results as in Sec. IV B. The calculation at finite p implies a shift in the zero-momenta of the propagators and 
derivative vertices, p° ^ p° ± /x [28, 30], depending on the charge states of the particles. The correction to log Z{p) 
from diagram (a) in Fig. 5 with the interaction from Eq. (12) is given by 



logZ(,)(/.) = — ^ f3V 3 {V+ ^V-Y + mi {U+ + C/_) {AD + [/+ + [/_) (26) 
^ ' 16 L J 



2 



\2 



(b) 



(c) 



(d) 



FIG. 5: Correction to logZ(/i). Diagram (a) shows the tttt interaction in the heavy p limit, diagrams (b)-(d) the interaction 
via explicit vector meson from Eqs. (8) and (9). 
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FIG. 6: Corrections to rrid or CF. Dashed-dotted line: result from the gauge invariant subset of diagrams from the first two 
columns of Fig. 3. Dashed line: Result from diagrams (b) + (c) from Fig. 5. Dotted line: heavy p limit from Sec. IV B. Solid 
lines: free tt gas, free p gas, and the ppmr interaction from Fig. 5 (d). 



with 

oo ^ oo 

U± = X f dk — n\w ± /i], V± = 4t / dk k^ nVw ± d (27) 
J uj J 



and D from Eq. (17). Applying Eq. (18) to logZ(a)(/J.) reproduces the result for the photon selfenergy in the heavy 
p limit from Eq. (24) which is shown to be gauge invariant in Appendix C 4. 

Thus having established that equivalence of photon selfenergy and charge fluctuations (Eq. (3)) holds on the 
pcrturbative level, we are encouraged to evaluate the diagrams of Fig. 4 by differentiating the appropriate terms in 
log Z with respect to the chemical potential. The diagrams for log Z corresponding to the photon self energies given 
in Figs. 3 and 4 arc displayed in Fig. 5(b,c,d). (Details can be found in Appendix C3). 

In Fig. 6 corrections to the electric mass of a free pion gas due to different sets of diagrams are shown. As a 
reference, we also plot the results for gases of noninteracting pious and noninteracting p- mesons (" free tt" and " free 
p"). The electric mass from the diagrams of Fig. 2 with the tttt interaction in the heavy p limit is plotted as the 
dotted line. The electric mass from the diagrams in the first two columns of Fig. 3 with dynamical p is plotted as the 
dashed-dotted line. At low temperatures, both results coincide (in detail this is also plotted in Fig. 16). However, at 
higher temperatures we observe significant differences which shows, thus, that the p obtains importance as an explicit 
degree of freedom. 

The diagrams (b) and (c) from Fig. 5 correspond to the first two columns of Fig. 3. Additionally, they provide 
photon selfenergies with ^pp and ^"fpp vertices from Fig. 4, diagrams (2a), (3a), (7a), and (8a). As shown in Fig. 6 
(dashed line), these additional 7/9 couplings obtain some minor influence above T ^ 150 MeV. 

Additionally, in Fig. 4 there are diagrams with ppmr couplings from Eq. (9). The diagrams (lb), (2b), (3b), (4b), 
and (7b) correspond to diagram (d) in Fig. 5. In the heavy p limit these diagrams do not contribute. However, for 
dynamical p-mesons these diagrams contribute significantly due to the sum over the spin of the p. In Fig. 6 the 
resulting electric mass is displayed as the solid line (" ppTTTr" ) . 
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VI. RELATIVISTIC VIRIAL EXPANSION 



In Ref. [28] the electric mass has been determined using chiral tttt interaction and thermal loops leading to results 
that show large discrepancies to a virial calculation of m^j. Before we discuss these differences in Sec. VIA, VI B let 
us review the theoretical framework first. The virial expansion is an expansion of thermodynamic quantities in powers 
of the fugacities e^'', while the interaction enters as experimentally measured phase-shifts. Consequently, all orders of 
the interaction are taken into account. Thermal loops, on the other hand, respect quantum statistics (Bose-Einstein 
in our case) and, thus, contain an infinite subclass of the virial expansion. However, the interaction only enters up 
to a given order. Thus, the loop and virial expansion represent quite different approximations and it will depend 
on the problem at hand which is the more appropriate one. The effect on quantum statistics can be considerable. 
For example at T = 170 MeV the values for the electric mass m^j of the free tt gas or the two-loop diagrams, Eq. 
(24), change by 20% and 38% (!) respectively, if we take the Boltzmann limit. Therefore, it is desirable to have a 
density expansion that respects particle statistics as well as sums all orders of the interaction. While this might be 
very difficult if not impossible to do in general, it can be done up to second order in the (Bose-Einstein) density. 

The partition function can separated into a free and an interacting part. 



logZ = logZo + 4' 4' b{h,i2) 



(28) 



in an expansion in terms of the chemical potential fj, with 



exp{Pfij) for j — 1,2 the fugacities. In the 5*- 



matrix formulation of statistical mechanics from Ref. [42] the second virial coefRcient 6(ii, 12) can be separated into 
a statistical part and a kinematic part containing the vacuum S'-matrix according to 



b{ii,i2) 



V 
47ri 



(2^ 



dE e-'^^vTOF T 



AS-\E)^S{E) 



(29) 



where A is the (anti)symmetrization operator for interacting (fermions) bosons and the trace is over the sum of 
connected diagrams (index "c"). In Eq. (29), V is the Volume, k is the momentum of the n-particle cluster in the gas 
rest frame and E ~ s^l"^ stands for the total cm. energy. The labels ^1,12 indicate a channel of the 5-matrix with 
i\ -\- 12 particles in the initial state. For the second virial coefficient, ii = 12 = 1. 

For TTTT scattering, Eq. (29) can be integrated over k and the 5-matrix can be expressed via phase shifts, weighted 
with their degeneracy [28]. With B2 = b{ii,i2)/V in the limit V 00 one obtains 
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(30) 



where the second line has been obtained after integration by parts (assuming (5^ — s- as E' — > 2m.,^). The sum over 
phase shifts 6^ (isospin /, angular momentum t) is restricted to 1 + 1 — even and Ki are the modified Bessel functions 
of the second kind. The virial expansion in this or similar form has been applied in numerous studies of the thermal 
properties of interacting hadrons as, e.g., [43, 44], among them the electric mass [28]. Note that, e.g. in [28], Bose- 
Einstein statistics is taken into account for the non-interacting, free gas, part. However, it is also possible to include 
particle statistics for the interacting part. This means the summation of so-called exchange diagrams as outlined in 
Ref. [42], Sec. VIIB. We employ this idea and also include a finite chemical potential. This is achieved by projecting 
the binary collisions of pions in different charge states to the isospin channels [28]. Additionally, the interaction T 
matrix is boosted from the gas rest frame to the two-particle cm. frame [65] and the T-matrix is defined via phase 
shifts with the final result 



tttt), Bosc 



= 4^ 



dE dx dk 



Ek^ 



X/S2 + fc2 



5l{E) (n[cJi + ^i\n[ijj2 + /i] + n[bJi - /-i]n[cj2 - yu]) 



2m„ -1 

-I- 6l{E) {n[u}i + ^]n[w2] + n[wi — /i]n[w2]) + 3 6\{E) (n[wi -I- ^]n[a;2] 4- n[uji — /i]n[ci;2]) 
-f 5l{E) ^^"-['^i + ^-l\n[uJ2 - l-i] + ^'n'[uji]n[uj2]^ + 3 Sl{E) n[LJi + fi]n[uj2 - ^] 

/2 1 

+ ^oi^) ( -^n[uji + fj]n[uj2 - fJ.] + -n[uji]n[uj2] 



(31) 
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A more explicit derivation of this result can be found in Appendix D. The first line of Eq. (31) corresponds to tttt 
scattering with a net charge of the tttt pair of |C| = 2, the second line to |C| = 1 and the third and 4th line to C = 0. 
The boosted Bose-Einstein factors which arise after summations over exchange diagrams arc 

1 ^1 ^ , kQx \ (I kQx 



k 



2 



^f-['-WTl^) ' Q=2V/^"^ (32) 

with the momentum of the pion Q = Qc.m. in the two-pion cm. frame. 

Obviously, the chemical potential can not be factorized in Eq. (31) so that the expansion is rather in powers of 
Bose-Einstein factors n than in powers of e^^ as in a conventional virial expansion. Eq. (31) contributes also to 
higher virial coefficients. The situation resembles the case of a free Bose-Einstein gas that contributes to all virial 
coefficients which can be seen by expanding the Bose-Einstein factor in powers of e^^ . Therefore, in the following we 
will refer to the expansion (31) as "(low) density expansion". The term "virial expansion" will be reserved for the well 
known expansion in terms of classical (Boltzmann) distributions. We note, that in the Boltzmann limit the standard 
expression for the virial coefficient, e.g. Eq. (9) of Ref. [28] is recovered; setting additionally = we obtain Eq. 
(30). 

The connection of B2 (p-) to physics is given by 

\ogZ{fi) = VB2{fi), P{f,)^^, ,nl^e'(^] (33) 



13 ' \dfj: 



ti=0 



where P is the correction to the pressure Note that for the electric mass the contribution ^ 6q vanishes in the 
Boltzmann limit (and is small anyways). The form of Eq. (31) makes it as easy to use as the common virial 
expansion, inserting the tttt phase shifts Sq, Si, and 5q which we adopt from Ref. [44]. The inelasticities of the tttt 
amplitude arc small in the relevant energy region and we have not taken them into account in Eq. (31). 



A. Density expansion versus thermal loops 

It is instructive to see to which extent the thermal loop expansion and the extension of the virial expansion from 
Eq. (31) agree. To this end we need to match both approaches by extracting the scattering amplitude from our model 
Lagrangian and insert it into Eq. (31). For simplicity, we first study the tttt interaction in the heavy p limit at 
and evaluate Eq. (31). As this interaction is not unitary, one has to go back to the original ^-matrix formulation and 
express it in terms of the (on-shell) T-matrix [42] which can then be calculated from theory. Given the normalization 
of the T-matrix used in this paper, 5 = 1— T, the right hand side of Eq. (29) can be written as 



dE dE ) SttOE 



^{T + T^) 



+ (9.T] (34) 

647r2 \e dE \E ' ^ ' 



Using the relation between ^-matrix and phase shifts, S ~ e^'*, we find 

where the connection between isospin amplitudes and their projection into partial waves T/ is given in Eq. (B5). 
Inserting this expression into Eq. (31) leads to the density expansion based on a given model amplitude. We note 
that the second term in Eq. (35) is quadratic in the amplitude and vanishes for real amplitudes. Therefore, close to 
threshold, where the amplitudes arc small and real, the quadratic term can be neglected. However, with increasing 
energy unitarity requires that the imaginary part of the amplitude will become sizable so that the second term cannot 
any longer be neglected. This is especially the case if the amplitude is resonant. Consequently, the use of point- 
like interactions at tree level which are always real and not unitary might lead to rather unreliable predictions for 
thermodynamic quantities. 

Before we discuss the importance of unitarity, let us first establish that the density expansion of Eq. (31) and the 
loop expansion lead to the same results if both methods are based on the same point-like interaction. The partial 
amplitudes T/ for the tttt interaction in the heavy p limit are obtained from Eq. (B4) by neglecting s, t, u, and T in 
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FIG. 7: Electric mass from dynamical p exchange. Solid line: from the diagrams in the first two columns of Fig. 3. Dotted line: 
Bose-Einstein density expansion from dynamical p exchange (no imaginary parts, Fp — > 0). Dashed line: Same, but Fp = 150 
MeV. 

the denominators and s = E'^. Inserting the result in Eq. (31) and calculating the pressure from Eq. (33) we obtain 
exactly the same result as for the thermal loops from Eq. (26) at /i = 0. We have also verified that this agreement 
holds in a simple 0^ theory of uncharged interacting bosons. Calculating the electric mass in both approaches for the 
TTTT interaction in the heavy p limit (Eqs. (24,3) and (31,33)), we again find perfect agreement. 

Consequently, and not so surprisingly, both thermal loop and density expansion lead to the same result, if the 
interaction in the density expansion is truncated at the appropriate (unitarity violating) level. This is also true in 
the classical (Boltzmann) limit. In this limit, a similar equivalence has been found in [28] using an effective range 
expansion for the amplitude; see also [45] for a related equivalence for propagators. 

While it is comforting to see that both approaches agree in the same order of density and interaction, this agreement 
highlights a possible problem for the loop expansion. If the order of the interaction considered violates unitarity the 
second term of Eq. (35) is ignored and the loop expansion may lead to unreliable results for the pressure etc. This is 
of particular importance if the amplitudes are resonant, as it is the case for the p-exchange. 

In order to see these effects we concentrate on the gauge invariant set of diagrams given in the first two columns 
of Fig. 3. The result for these diagrams is given in Eqs. (C2, C4) and plotted in Fig. 7 as the solid line. In the 
calculation of these thermal loops we have made the following approximations, see Appendix CI: (I) The poles of 
the p have been neglected in the contour integration (see the explanation following Eq. (C8)). (II) The p has no 
width, i.e. the p propagator is given by Z?'^" from Eq. (22). (Ill) Only the real parts of the thermal loops have been 
considered. 

In the following we test these approximations by comparing the thermal loop result with a suitable "toy model" 
low density expansion. For the interaction driving the low density expansion we take the partial waves from Eq. (B4) 
and project out the T/ by the use of Eq. (B5). Furthermore, we set Fp = in Eq. (B4) in this interaction. Third, 
we consider only the term linear in T in Eq. (35) for the density expansion. This means that imaginary parts are 
neglected. The low density expansion, constructed in this way, exhibits the same approximations (II) and (III) as the 
calculation of the thermal loops above, i.e. the zero width and the reduction to the real part only. The result of this 
"toy model" low density expansion is plotted in Fig. 7 as the dotted line. 

Both the results from thermal loops (solid line) and the density expansion (dotted line) agree closely. The small 
deviation of both curves is due to the additional approximation (I) which we have made in the calculation of the 
thermal loops, i.e. neglecting the poles in the contour integration. Note also that other partial waves than Tq, T^, 
and are present in the results from the thermal loops because the p exchange contains all partial waves. However, 
from the agreement found here, we may conclude that these higher partial waves give negligible contributions (at least 
in the present p-exchange model). 

In our "toy model" low density expansion, we can allow for a finite width in the p-propagator. This implies that the 
p-propagator is given by Dp = [p'^ -m'^p + impT{y/s)]~^ , where T{y/s) = F(TOp)(mp/£^^)(i?^ - 4mJ)^/^/(rnp - 4mJ)^/^. 
With this modification, we evaluate again the electric mass. However, as we still use only the term linear in T in Eq. 
(35), any imaginary parts of the amplitude arising from the finite width are still ignored. The result is shown as the 
dashed line in Fig. 7; the electric mass hardly changes. 

Let us now discuss the effect of the imaginary parts of the amplitude. To simplify the discussion let us restrict 
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FIG. 8: Main plot: Contribution to rn^i in the = (1, l)-channel of tttt scattering. Solid line: From experimental phase 

shift. Dotted lines: Terms linear and quadratic in T from Eq. (36) with T from a unitary p exchange model. Dashed line: Sum 
of linear and quadratic term. Gray line: free p-gas (Boltzmann). 

Insert: With T from the LO chiral Lagrangian. Dashed line: Tree level. Solid line: Unitarization with A"-matrix. 



ourselves to the vector-isovector (/, J) — (1, 1) channel, which is dominated by the p-resonance. We will also work 
in the Boltzmann limit as effects due to unitarity are independent of the statistical ensemble. The model amplitude 
is simply the s-channcl p-cxchangc diagram with a p-propagator as given above. This amplitude is unitary by 
construction and describes the scattering data in the {I, J) = (1, l)-channel well (see Fig. 15). With the (complex) 
T-matrix the electric mass TOqI is given by 

OO OO t-> 

ml{n = 0,I = e=l) = -f- / dE Q E KiiPE) ReTl + / dE K2{(3E) {ReT^ ImT^) 

(36) 

where the first term is linear and the second quadratic in the amplitude. It is the second, quadratic, term where the 
imaginary part of the amplitude enters. In Fig. 8 the different contributions to the electric mass according to the 
decomposition Eq. (36) are plotted. As a reference we also show the result using experimentally measured phase shift 
6\ (solid black line). Obviously the contribution from the quadratic term {^^Tl, quad.") is dominant, and adding the 
linear ("T^^, lin.") and quadratic terms wc obtain good agreement with the result from the 5\ experimental phase 
shift. This is to be expected as Tl fits vacuum data well. Note that the linear term alone vastly undcrpredicts the 
electric mass. Thus the imaginary part of the amplitude is essential for the proper description of the fluctuations. 

Furthermore, the electric mass from the free p-gas (Boltzmann statistics) (gray line. Fig. 8) agrees well with the 
results from the experimental phase shift and unitary p-model. Indeed, it can be shown that the TTTr-interaction via 
unitary s-channel p-exchange in the limit of vanishing width leads to a contribution to log Z equal to that of a free 
p-gas [44, 46]. In this limit 5\{E) = ttQ{E — nip) allowing for an explicit evaluation of Eq. (31) in the Boltzmann 
approximation. For Bose-Einstcin statistics the situation is more complicated. In Ref. [46] it has been shown for 
mcson-baryon interaction that also in this case the interaction of two particles via a narrow resonance TV* leads to the 
same grand canonical potential as from a free 7V*-gas with the corresponding Fermi-statistics for A^*; however, the 
proof requires a self-consistent medium modification of the N* width and a consideration of larger classes of diagrams. 

While in our toy model we could simply restore unitarity by introducing a p width, in a more complete calculation 
this is considerably more difficult. For example, using a p-propagator with a width in the diagrams of Figs. 3 and 
4 leads to additional photon couplings to the intermediate pion loops, which generate the p-width. This is simply 
a consequence of gauge invariance (see Appendix C4). Therefore, introducing unitary amplitudes while maintaining 
gauge invariance is a non-trivial task. 

An alternative approach to assess the role of unitarity is to unitarize a given amplitude using the A'-matrix approach 
(see, e.g. [47]). This approach docs not add any additional dynamics, and therefore provides a good estimator on the 
importance of unitarity alone. Using the AT-matrix approach we can in principle take any of the interactions discussed 
in this paper. Here we choose the interaction in the (1, l)-channel from the LO chiral Lagrangian given in Eq. (6). 
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FIG. 9: CF or electric mass for the interacting tvtt system. The £S, /li^J calculations are from Ref. [28] and the conventional 
virial expansion (Boltzmann statistics) reproduces results from Ref. [28]. 

Details of the calculation can be found in Appendix B2. Maintaining gauge invariance in a A'-matrix unitarization 
scheme requires special care and is beyond the scope of this paper. Ignoring this issue, we can compare the electric 
mass from the unitarizcd version using Eqs. (B9,31,33) with the tree level amplitude using from Eq. (B8) and then 
Eqs. (35,31,33). The results arc plotted in the insert of Fig. 8 and show only a small correction due to unitarization. 

Consequently, unitarity by itself is not as crucial as the dynamics which generates the resonance. In other words 
as long as the phase shift is slowly varying with energy unitarity corrections are small. A resonant amplitude on the 
other hand corresponds to a very rapidly varying phase-shift. Since it is the derivative of the phase-shift which enters 
the density expansion, resonant amplitudes are expected to dominate. Consequently, a resonance gas should provide 
a good leading order description of the thermodynamics of a strongly interacting system. 

Note that the unitarized amplitude T^^^ ^ from Eq. (B8) corresponds to a unitarization via the Bethe-Salpeter 
equation in the limit where the real parts of the intermediate TTTr-loops arc neglected; the freedom in the choice of the 
real part (loop rcgularization) can be used to fit to experimental phase shifts, which in turn introduces the missing 
dynamics (see, e.g., Rcfs. [48, 49]) This should lead to more reliable predictions [50]. 

To conclude this analysis of the density expansion, it appears that the low density expansion of Eq. (31), using 
experimental phase shifts, will give the most reliable results, while a simple hadron gas calculation should provide a 
reasonable first estimate for the fluctuations of a system. Finally, there are certain features of the p model from Sec. 
V which can not be taken into account in the low density expansion: The ppmr and {'j)^pp interactions discussed in 
Sec. V (Fig. 4) arc a consequence of the p being treated as a heavy gauge particle; these features will be missed in 
the low density or virial expansion in which the p is not more than a resonant structure in the tttt amplitude. These 
considerations will be taken into account in the final numerical result from Sec. VIII. 



B. Numerical results for the interacting pion gas 

In Fig. 9 the results so far obtained are compared to Ref. [28] (gray dashed lines). The electric mass for pions 
interacting in the heavy p limit from Sec. IV B is indicated with the dotted line. Taking into account that the 
interaction from Eq. (12) is around 3/2 times stronger than the one from the LO chiral Lagrangian, the calculation is 

(2) 

consistent with the Ct„t calculation from Ref. [28] which we have also checked analytically. The result for dynamical 
p exchange (black dashed line) contains the contributions from free pion and p gas and the diagrams from Fig. 5 (b), 
(c), and (d). The difference to the heavy p limit shows the importance of the p as an explicit degree of freedom in 
the heatbath. 

Up to r ^ 130 MeV the dynamical p exchange contributes with the same sign as the virial expansion from Ref. 
[28] although they differ largely in size due to the lack of imaginary part in the loop calculation, especially in the 
(J, £) = (1, l)-channcl as discussed above. Also the p-modcl docs not describe the = (2,0) amplitude very well. 

For the low density expansion from Eq. (31) and the virial expansion from Eq. (30) we use the phase shifts from 
Ref. [44]. Note that there is a partial cancellation from the Si and Sq partial waves [28]. 

The calculation from Ref. [28], which of course contains also the contribution, shows a very distinct 
result. The reason is twofold: on one hand, unitarity is not preserved (see discussion in Sec. VIA). On the other 
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FIG. 10: Ds from Eq. (4) for the interacting tttv system. The heavy p hmit from Eqs. (24,25) includes the contribution from 
the free tt gas. The contributions from free tt and free p gases are included for the result from Eqs. (C9,C11,C12), the latter 
indicated as "dynamical p". The virial and density expansions are indicated according to their underlying particle statistics 
with "Boltzmann" and " Bose-Einstein" , respectively. 

hand, the thermal loops in the C^J calculation pick up high cm. momenta where the theory is no longer valid and 
the dependence of the NLO interaction on high powers of momenta introduces artifacts. Note that the size of the 
correction from C^^^ alone is larger than the one from C^} for T > 80 Mev. 

The results for the observable Ds from Eq. (4) are displayed in Fig. 10. Corrections to the entropy are included: 
from Eq. (25) for the heavy p limit, from Eq. (C2) for the case with dynamical p and from Eq. (33) for the low 
density expansions. For the thermal loops, indicated by "models with p exchange", Ds is suppressed. This is due to 
the large negative correction to {SQ^') as has been seen in Fig. 9. The virial expansion and the density expansion 
coincide better with each other than in Fig. 9 and can be roughly approximated by a gas of noninteracting pions and 
rhos. 

Having contrasted virial expansions and dynamic p model in Sec. VIA, the most realistic results for CF and Ds 
for the interacting tttt system are given by the Bose-Einstein density expansion from Eq. (31). While this result is 
somewhat below the estimate of a gas of free pions and p-meson, it is nowhere near the value of Z?^'^^ ~ 0.034 for 
the quark gluon plasma. 

VII. HIGHER ORDER CORRECTIONS 

Both the density expansions and p models from the last sections are quadratic in density, i.e., the statistical factor 
n. However, at the temperatures of the hadronic phase higher effects in density play an important role. Virial 
expansions become complicated beyond the second virial coefficient and no experimental information exists on three 
body correlations. Performing resummations is, therefore, of interest. This will include the density and strong coupling 
g to all orders. Of course, this can not be done in a systematic way; resummations only contain certain classes of 
diagrams at a given order in perturbation theory. In all resummations, log Z is calculated at finite p and then Eq. 
(18) is applied in order to obtain the electric mass. We have convinced ourselves in Sec. V that this is a charge 
conserving procedure. 

We start with two natural extensions of the basic interaction diagram (a) in Fig. 5, displayed in Fig. 11 (n) and (r) 
using for both of them the effective interaction of the heavy p limit from Eq. (12). Alternatively, one can use the LO 
chiral Lagrangian from Eq. (6). As found in Sec. HI, results for the dominant part of this interaction are obtained 
by simply multiplying g in the following by a factor of 2/3. However, one should keep in mind the unitarity problems 
of these simplified point-like interactions which have been addressed in Sec. VIA. 

For the calculation of diagram (n) we utilize an equation of the Faddeev type. The Faddeev equations, usually 
used in three-body scattering processes as in Ref. [49] in a different context, are an easy way to sum processes whose 
elementary building blocks are of different types, as in this case loops of neutral pions with chemical potential p = 
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FIG. 11: Resummation schemes: necklace (n) and ring (r). Below, the tadpole medium correction of the p propagator is 
displayed (t). 



and charged loops: 

log ^(n) (m) = ^ Q ^± + ^± + ^o) 

b± = g'c±+ g'l±{b±+bo) 

with g' = —g^/m^p. The first loop in the chain is labeled a, the last one c, and I means an intermediate loop. The 
indices "±" and "0" label charged and uncharged loops, respectively. It is instructive to expand Eq. (37) loop by 
loop which shows that the structure indeed reproduces all sequences of charged and uncharged loops, of all lengths. 
There is a symmetry factor of 1/2 for every loop of neutral pions and a global factor of 1/2 for every pion chain. The 
solution of Eq. (37) is found in Appendix E. The result of resummation (n) is plotted in Fig. 12 together with its 
expansion up to g^/m^ (dashed line) and up to g^/m^ (dotted line). 

The summation (r) of Fig. 11 with the interaction from Eq. (12) exhibits a symmetry factor of 1/A^ for a ring 
with TV "small" loops (see Fig. 11) which after summing over leads to the occurrence of a logarithmic cut in the 
zero-component p° of the momentum of the "big" loop. Due to this obstacle for the contour integration method [30], 
usually only the static mode p*^ = is calculated, although new studies overcome this problem [51]. In the present 
approach, we can calculate the ring with N "small" loops explicitly before summing over N. This avoids, thus, the 
problem of the logarithm at the cost of having to cut the series at some Ny^ax- On the positive side, all modes are 
included, and not only the p° = static contribution. The result up to eight "small" loops has already converged up 
to r ~ 200 MeV and is displayed in Fig. 12 as (r). The explicit solution can be found in Appendix E. 

Note that in the resummation schemes we do not consider the vacuum parts of the loops, i.e. we do not renormalize 
the vacuum amplitude. This excludes potential double counting issues in the final numerical results in Sec. VIII 
where resummations and density expansion are added: Renormalizations of the vacuum amplitude are supposed to 
be be included in the phase shifts that are used in the density expansion. 

There is an additional resummation scheme that sums up the ppmr interaction required by the gauge invariance of 
the p (sec Eq. (9)): One can consider diagram (b) and (c) of Fig. 5, dress the p propagator as indicated in (t) of Fig. 
11, and finally take the heavy p limit as in Sec. III. This leads to the same result as a renormalization of the static p 
propagator — l/m^ of diagram (a) in Fig. 5 for the tttt interaction in the heavy p limit: The resummed pion tadpoles 
can be incorporated by a mass shift, 

2 2 

ml±-.ml + ^{U+ + U^+2D), m^,, ^ + |-([/+ + [/_) (38) 
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FIG. 12: Left panel, main plot: Resummations (n) and (r) from Fig. 11 and their expansions up to /m^ (dashed hne) and up 
to /m^p (dotted hne) which are the same for (n) and (r). Insert: resummation (t), for order g^/m"^ and higher. Right panel: 
CF over entropy, Ds- Resuh for heavy p limit and dynamical p as in Fig. 10. The result including the resummations (see text) 
is shown as the dashed dotted hne (S'). 



for charged and neutral p. The contribution to md from this modification is shown in the insert of Fig. 12 as (t). 
The thermal po mass from Eq. (38) at yu = is = 824 MeV at T = 170 MeV which is slightly more than in other 
studies [52] . This is certainly due to the omission of the p — > tttt —f p selfenergy which also contributes and is required 
by the gauge invariance of the p-meson. In the counting of the present study, the p tttt p selfenergy is statically 
included in the resummation (n) of Fig. 1 1 . 

To the right in Fig. 12 the normalized CF over entropy, Ds from Eq. (4), are plotted. For comparison, the result 
at from the dynamical p-exchange (see Fig. 10) is shown with the dashed line. We include now the resummation 
(n) but only with three or more loops, or in other words, at g'^ and higher in the interaction in order to avoid double 
counting with the g^ contribution. We have already seen in Fig. 12, left panel, that both resummations (n) and (r) 
contain the same diagram at order (linear chain of three loops). Thus, again in order to avoid double counting, 
we include the resummation (r) requiring at least three of the "small" loops, see Fig. 11; this means that only 
contributions of order g^ and higher are included. Finally, we add the resummation (t) including the orders g* and 
higher, which again avoids double counting of the g^-contribution. Summing in this way the resummations to the 
g^-result (dashed line) for both ISQ"^) and S, the resulting Ds ~ (SQ^) / S is indicated as S' with the dashed-dotted 
line in Fig. 12. 

The resummations have a large effect on ((5Q^) (see Fig. 12, left) whereas their effect on the entropy is much 
smaller; the entropy is efficiently suppressed for higher orders in the coupling. This explains, why the result E' shows 
such a large difference compared to the results at order g^ (dashed line). 

For the resummations (n) and (r), we have ensured that we recover the results from Eqs. (24) and (26) at the same 
order of the interaction. We have also verified that the results from Ref. [53] at external momentum p of the p being 
zero (p" = 0, p ^ 0) match the p self energies at /_i = that are implicitly or explicitly contained in the resummations 
(n) and (t). 

A possible extension of the diagrams discussed here is given by resummations of super-daisy type: the pion propa- 
gator is dressed by a series of pion tadpoles; the propagator of the tadpole loop itself is again dressed which constitutes 
a self consistency condition. E.g.. this leads to a thermal mass of the pion ra-jr ^ 170 MeV at T ~ 170 Mev. However, 
one should realize that the lower orders in the coupling g of a super-daisy expansion are already covered by the 
resummations considered before: it is easy to see that the super-daisy resummation introduces additional diagrams 
only at order g^ and higher {g^ and higher for resummation (t)) and, thus, can be neglected. 



A. Extension to SU{'A) 

In order to obtain a more realistic model for the grand canonical partition function, the leading contributions from 
the interaction of the full SU (3) meson and vector meson octets is considered. Obviously, the leading contribution 
to the CF from strange degrees of freedom is simply the free kaon gas. Here we want to discuss corrections due 
to interactions of kaons with pions. The most important of those is the resonant p-wave interaction involving an 
intermediate AT* (892) meson. This is quite analogous to the p meson in the tttt case, discussed previously. The $ 
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FIG. 13: Corrections to the electric mass or CF, (SQ^) /{e^VT^) for yriC interaction. The density and virial expansions are 
from Eqs. (F4) and (F6), respectively. The loop expansions "ttT^ dynamical" and "nK contact" are from Eqs. (F1,F2) and 
(F3), respectively. The solid line shows the electric mass of a gas of free ft;(800), -ft'* (892), _K'o(1430), and _ft'J(1430) mesons. 

meson, on the other hand, only enters if interactions between kaons are considered. These are sub-leading as pions 
are more abundant and thus -kK interactions are more important. 

As in the previous sections we describe the meson-meson interaction by dynamical vector meson exchange, second, 
by an effective interaction, and, third, by rcaUstic phase shifts via a rclativistic Bosc-Einstcin density expansion. For 
processes which contain at least one pion, the dynamical vector meson exchange is mediated by the A'* (892). The 

effective contact interaction is taken from the LO chiral meson-meson Lagrangian in its SU (3) version, for the 
Trif -interaction. The density expansion of nK scattering is obtained following the same steps as in Sec. VI. Details 
of the calculations are summarized in Appendix F. 

In Fig. 13 the CF from the different models are shown. For the virial and density expansion the phase shifts 

1/2 1/2 1/2 

have been taken from the parametrization of Ref. [54] for the attractive channels 5^ ,5i ,82 , corrected for the 
parameters of the iirg(1350) resonance (nowadays, ^0(1430) in the PDG [55]) as reported in Ref. [43]. The repulsive 

phase shift is from Ref. [56]. The phase shifts plotted in Fig. 4 of [43] up to s^/^ = 1 GeV have been reproduced. 

The situation resembles the case of tttt scattering from Fig. 9: Thermal loops with dynamical vector exchange 

(2) 

or with effective interaction via >C^j^ show large discrepancies to the virial and density expansions, this time even 
more than in the tttt case. The reasons are similar as those found in Sec. VIA: The repulsive {I,t) = (3/2,0) 
partial wave is not well described by ttK scattering via AT* (892) and unitarity problems of the thermal loops show 
up. The contributions from both the virial expansion and the low density expansion are large compared to the virial 
corrections in the tttt sector (see Fig. 9). This seems surprising as in the ttK system the kaon has a large mass which 
should suppress contributions kinematically. However, in the considered channels of i:K scattering, four resonances 
are present, k(800), /v*(892), A'q(1430), and A'2(1430) and we know from Sec. VIA that resonances give a large 
positive contribution to mci[66]. The electric masses from these resonances, treated as free gases (Boltzmann), is 
plotted in Fig. 13 with the solid line. We find the same pattern as in the discussion of Fig. 8 for the free p: the 
virial corrections from resonant phase shifts are well described by a free gas of the same resonances. Furthermore, the 
repulsive (5q^^ phase shift is very small. 

As the outcome for the density expansion in Fig. 13 shows, the inclusion of Bose-Einstein statistics is important 
(compare to the virial expansion which uses Boltzmann statistics) . We consider the density expansion to provide the 
most reliable prediction. 

VIII. NUMERICAL RESULTS 

In the discussions in Sees. VI A and VII A good reasons have been found that at quadratic order in density n the 
Bosc-Einstcin density expansion gives the most realistic results. For the final numerical results we include therefore 
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FIG. 14: Final results for charge fluctuations (electric mass) and Ds- Results of the Bose-Einstein density expansions with the 
dashed-dotted lines. Adding pp-K-K and K* K'-ktt contributions, the resummations, and free mesons up to 1.6 GeV, the results 
are indicated with the dashed lines. For comparison, m^i and Ds from free mesons alone, without any interactions, are also 
plotted (dotted lines). 



the TTTT and the ttK density expansion from Eqs. (31) and (F4). The dashed-dotted lines in Fig. 14 show electric mass 
and normalized charge fluctuations Ds from Eq. (4) for the sum of the two density expansions. At order n^, there 
are additional photon selfeirergy diagrams with (7)7^^ and ppin: vertices from Fig. 4. As discussed at the end of Sec. 
VI A these contributions are not included in the density expansion but a consequence of the p being introduced as a 
heavy gauge field. The same applies to the K*K*Tnr diagram (d) from Eq. (F2). Thus, we include these additional 
contributions for ((5Q^) and S. 

At higher orders in density one has to rely on resummation schemes. Including the resummations in the final results 
does not to lead to double counting: Resummations at order g'^ and upwards in the strong coupling correspond to 
diagrams with three and more loops and, thus, to contributions higher than quadratic in density. We include (with 
5^ and higher) the summations (n), (r), and (t) from Sec. VII. Note that for TOei there is a partial cancellation of 
sizable contributions from the resummations and the {'y)jpp, ppmr, K*K*'!nT diagrams. 

In order to obtain a more realistic picture, we also include as free gases all mesons from the PDG [55] which have 
not been considered so far, up to a mass of 1.6 GeV. Note that we do not add free mesons that have the same quantum 
numbers as the density expansions namely cr(600), /9(770), k(800), K*{892), K^{U30), and /q(1430). We have seen 
in Sec. VIA that their coirtribution via phase shifts in the density expansions is roughly of the size as if they had 
been included as free particles. Adding all coirtributions mentioned, the results are iirdicated with the dashed lines in 
Fig. 14. 

Compared to the density expansions the final results do not change much. The influence of heavier mesons than 
those considered in this study is, thus, well controlled. Many of the heavier resonances that have been added here 
as free gases are axials which decay into three particles. To include them in a density expansion would require the 
consistent treatment of three body correlations. 

Concluding, we can assign Ds — 0.09 for temperatures 120 < T < 200 MeV which coincides (incidentally) quite 
well with the result if one simply considers free, noninter acting, mesons up to masses of 1.6 GeV. The latter case is 
indicated with the dotted lines in Fig. 14. 

Theoretical uncertainties in the present study arise from the omissioir of diagrams such as the (small) eye shaped 
diagram mentioned in Ref. [28] already at g^. Furtherirrore, both resummations and density expansions are incomplete 
as they only partly include the in-medium renormalization of the resonances which drive the nresoir-mcson scattering, 
such as the (7(600), /o(980), or the p(770) itself [46, 53]. In this coirtext one can think of a more complete microscopical 
model: We have found in Sec. V and VI A that unitarity and a good description of the vacuum data up to high energies 
and in all partial waves are important. Such models exist, e.g., the chiral unitary approach from Ref. [48]. The 
medium implementation of such a model has been done in a different context, see [57, 58] and references therein. A 
generalization of the virial expansion from Ref. [59] to finite chemical potential and including Bose-Einstein statistics, 
as carried out here, would be feasible in principle. Such an ansatz [50] would allow to take simultaneously into account 
the medium renormalization of the (dynamically generated) resonances and the calculation of the grand canonical 
partition function at finite p as needed for a calculation of to^j. 
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IX. SUMMARY AND CONCLUSIONS 



For an estimate of charge fluctuations (CF) in the hadronic phase of heavy ion coUisions, we have calculated the 
effect of particle interactions. For the perturbative expansion up to two thermal loops, the tttt interaction has been 
described by vector meson exchange. The correlations induced by a dynamical p have been found significant by 
comparing to an effective theory where the p is frozen out. 

The photon self energies are charge conserving and shown to be equivalent to the loop expansion of the grand 
canonical partition function at finite chemical potential. We have pointed out that the inclusion of imaginary parts 
is essential for a proper description of the thermodynamics, especially if resonant amplitudes are involved. To second 
order in the density, it has been possible to include Bose-Einstein statistics in the conventional virial expansion. 
This "density expansion" can change the conventional results significantly. Moreover, for real amplitudes, we could 
show the equivalence of the loop expansion and the density expansion at all temperatures. However, the inclusion 
of unitary (complex) amplitudes is more straightforward in the density (virial) expansion. To the extend that two- 
particle correlations are dominant, the density expansion with Bose-Einstein statistics is, thus, the method of choice 
as it provides the same statistics as the thermal loop expansion and unitarity is automatically implemented by the 
use of realistic phase shifts. 

For an estimate of three- and higher particle correlations, a variety of summation schemes has been presented, 
all of which tend to soften the large first order correction of the thermal loop expansion. For the CF, higher order 
corrections have a large influence whereas higher orders for the entropy are small. 

For the CF over entropy, Ds, it has been shown that the influence of heavy particles beyond the interactions 
considered arc well under control; a final value of Ds — 0.09 has been found for temperatures 120 < T < 200 MeV. 
This result agrees quite well with the outcome from the free resonance gas, supporting the notion that resonant 
amplitudes dominate the thermodynamics. As lattice gauge calculations with realistic quark masses become available 
it would be interesting to see at which point these start to significantly deviate from a hadron gas. 
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APPENDIX A: FROM CHARGE FLUCTUATIONS TO PHOTON SELFENERGY IN SQED 

In this section, an outline for the proof of Eq. (3) for scalar QED is given. The argument follows Rcf. [30] where a 
similar connection is made for QED. If tttt contact interactions arc included according to Eq. (14), the steps outlined 
below are similar, but lengthier, and Ward identities for four-point functions have to be determined. 

CF are defined as (^Q^) = (Q^) ~ (Q)^! a-nd the expectation values are calculated via the statistical operator of 
the grand canonical ensemble with the charge chemical potential p = pq. One obtains immediately: 

{5Q^)=e'TV ^^{],) (Al) 

with Jq the zero-component of the conserved current, Q = / Jo = ^Jo- The expectation value of Jq — i(0*(9° -I- 
ieA^)(l) — (/)(9° — ieA^)(j)*) can be expressed in terms of the propagator 

fd{W\ ^ f d^ 



T E y (|;^2(p°-^)5(p°,p) (A2) 



\ dp J rp dp ^ J (27r)3 
where we have used p = eJ^ and the definition of the imaginary time propagator 

5„;3(xt;x't') = -Tt\pg T,[(/)A'a(xT)(/)t^^(x'T')]] (A3) 

where is the r-ordcrcd product in the modified Heisenberg picture, see, e.g., Ref. [60], and the Fourier transform 
is at equal times r, r+ and position x = x'. The /i-dependence of the propagator is given by 'p^ ~ iun ^ p where 
LUn = 27ri n T. With this, the derivative can be rewritten as 
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at zero chemical potential /i = 0. Using d/dp^Q = ~G{d / dp^Q~^)Q , the Ward identity in the differential form for 
scalar QED can be applied. The Ward identity connects the inverse propagator with the fully dressed vertex 
according to 

e^TT/(^^)^ = TV £ y'^[2e2t;(/,p)-e(2pO)c;(/,p)rO(p,c.„M/,p)" 

= rF(noo,„,t(fco = o,k^o) + n«'^,t(fco = o,k-.o)). (As) 



Factors of e and have been identified here with the bare ^^ixix and 77r7r vertices. In the step from Eq. (A4) to Eq. 
(A5) we have generated three propagators from one, and it should be noted that this takes place inside the momentum 
integral and summation. Therefore, the limit in Eq. (3) has to be taken before summation and integration. 

APPENDIX B: PION-PION INTERACTION 
1. Chiral nn interaction and vector exchange 

In this section the effective pion-pion contact interaction from Sec. Ill and its connection to the chiral Lagrangian 

(2) 

is discussed in more detail. For four pion fields the kinetic term of C)„t in Eq. (6) and the effective interaction in 
Eq. (12) have identical isospin and momentum structure. Comparing the overall coefficients leads to the result in Eq. 
(13) which differs from the KSFR relation by a factor of 3/2. Studying the low energy behavior of both theories helps 
solve this puzzle of the obvious violation of the phenomenologically well-fulfilled KSFR relation. The tttt amplitude 
at threshold from the LO chiral Lagrangian Eq. (6) and the effective interaction Eq. (12) is given by 



respectively which leads to the correct KSFR relation 

= (B2) 

This is due to the mass correction term proportional to A4 in Eq. (6). This term, however, does not have any 
momentum structure and immediately becomes small at finite pion momenta compared to the kinetic term. It has 
no inffuence in the results of this study. 

For finite pion momenta, higher order partial waves have to be included. We concentrate on the quantum numbers 
of the p-meson and obtain for tttt scattering via the LO chiral interaction in isospin / = 1: 

Tl^-^{t-u) (B3) 

which should be compared to the result from p-exchange from Eq. (8): 

s — u t — u t — s 



T^{dyn. p) = 5^ 



t — rn^ s — ni^ + impT{s) u — rn^ 



r2(dyn. p) = c/ ( + ±^ ] (B4) 



t — m?p u — rn^p 



where we have also given the result for for completeness, and is immediately obtained by crossing symmetry, 
T" = —2T'^. Projecting out the p-wave in both results (B3) and (B4) by using 

1 

Tl{s) = ^d{cose) P,(cos0) T'is,t,u) (B5) 
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FIG. 15: p-wave isovector tttt interaction. Dots; Partial wave analysis from Ref. [61]. Dashed line: cSl calculation. Dashed- 
dotted line: £1*^ calculation from Ref. [62]. Solid line: Effective interaction from Eq. (12). Thin solid lines: Explicit p exchange 
from Eq. (8) with and without (momentum dependent) width for the p. 



for (/, t) — (1, 1), making an expansion in p^ni, and comparing the coefficients, leads to the relation ■m'j^ — 4m^ — if^g^ 
which shows again the deviation of 3/2 from the KSFR relation up to a correction from the pion mass. However, taking 
only the s-channel vector exchange, which is given by the second term of Eq. (B4), we obtain after projection to 
the p-wave: 

ml ~ 4m2 = 2flg\ (B6) 

This is indeed the KSFR relation in Eq. (B2) with some small correction which vanishes when s is neglected against 
m?p in the denominator of Eq. (B4). Concluding, the restriction to s-channcl vector exchange in tttt scattering restores 
the KSFR relation in the p-wave expansion of the scattering amplitude. However, t- and u-channcl vector exchange 
is also present, and this leads to the effective interaction in Eq. (12) which is 3/2 times stronger than the interaction 
from the LO chiral Lagrangian. 

Fig. 15 illustrates the behavior of the different theories together with data from Ref. [61]: The LO chiral Lagrangian 

underpredicts the strength of the experimental amplitude. In contrast, the interaction up to C^^l and the effective 
interaction from Eq. (12) describe better the data at low energies. The explicit p exchange with width (thin line) 
delivers a good data description even beyond the p-mass. 

One more remark is appropriate in the framework of this section: In the treatment of the p-meson as a heavy gauge 
field, the covariant derivative introduces the up interaction as we have seen in Sec. HI. Additionally, the original tttt 
interaction from Eq. (6) remains in this derivative. In the present model, we have omitted this term, as has been 
also done, e.g., in Ref. [33]. This leads to better agreement with the data in the T^j^-channel and ensures the KSFR 
relation. It is possible to keep the original chiral interaction, but then additional refinements have to be added added 
as, e.g., in Ref. [37]. 



2. Unitarization of the TTTT-amplitude with the if -matrix 



The ii'-matrix is defined via the 5-matrix as 

Sk{E)^^^'^^\ K^^^ (B7) 
^ ' 1-iQK' 16ttE ^ ' 

with the tree level amplitude Ttrco from Eq. (B3) and the cm. momentum Q from Eq. (32). The unitarizcd amplitude 
^(«) 1 which is given by 

l + 2zOri/i?' 96^/2 ^^^^ 



can be parametrized via phase shift as 



1 -ReTi , 
5\ = - arctan (B9) 
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APPENDIX C: THE p-MESON IN THE HEATBATH 
1. Analytic results 



The analytical expressions and numerical contributions from the set of gauge invariant diagrams in Fig. 3 are given 
which are obtained from the interactions from Sec. III. With 



L± (a, b) := log 



vn? + (p ± qf — {b u! — ato'Y 



riijj + (p T 9)^ — {b Lo — auj')^ 



m'l + (p— qf — [buj + auj'Y 



m'^ + {p + qf - {boj + auj')^ 



(CI) 



where oj'^ = q^ + and cj'^ = + m^, we obtain for the real parts of the diagrams in Fig. 3, left column: 



n°i%)(fco = o,k^o) = 



e^g^ d d 



00 00 



{2^)ml da 9/3' 



|q=/3=i I dp I dq j"^ n[uj] n [a/ (3 J 
J J apujix)'-^ 



-8ml pq 



L_(a/Al)) 



n??a2)(fco = o,k^o) = 



e?g^ d d 



oc oc 



2(27rr da 9/3' 



a=/3=i I dp I dq 



pq 



n[l3ijj\ n [aijj'] 



(fco = o,k^o) 



X L+(a,/3), 



00 00 



n°5%(fco = o,k^o) = 



(27r) dc 



Bm^ pq 



6e^g^ d_ 



7^|q=i dp dq 



pq 



n\uj\ nlauj 







{2m^mpY - (m^ - {a^ - 1^^)^ L^{a, 1) 



00 OC 



a=i dp dq 



pq 

,,,,,'2 



nlujl n laoj 



X [-au {ml + (1 - a^)iu'^) L- {a, I) + u {ml + {a^ - 1 V^) L+ (a, 1)] 



n°e%(fco = 0,k^O) 



2^2 



oc oc 



-3e/g 



dp I dq n[uj] n [uj'] 



(2^)^m2 

X [- {ml -u? - iu'^) L- (1, 1) + 2ujLu'L+ (1, 1) 



logZ, 



(onl) 



3gV 

4T (27r)^ 



00 00 



dp I dq^ n[Lo] n [to'] [-8pq + {Ami - m^p) L- (1, 1))] 



(C2) 



In these expressions the poles from the p-propagator have been omitted as discussed in Appendix C2. The use of 
derivatives in Eq. (C2) cures infrared divergences which occur (sec Appendix C2). The logarithmic pole in the 
numerically relevant integration regions for p and q is in all cases is given by 



ml + ip-q f ~{bLJ + a uj'f = 



(C3) 



where a, b take values according to the arguments of L±{a, b) of Eqs. (C1,C2). The singularity leads to an imaginary 
part which wc neglect. The issue of imaginary parts is discussed in Sec. VIA. The diagrams from the second cohmm 
of Fig. 3 are calculated straightforward with the results 



ttOO 
ii(lcl) 



in the static limit (feg = 0, k — > 0). 



(5cl) 



'Hon. u^, 



7 ? 

— ^ 



(C4) 
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FIG. 16: Numerical results for the diagrams from Fig. 3, Eq. (C2) as a function of T [MeV]. Selfenergy Il/{e^T^) in [MeV"^] 
for all plots, except the correction to Z: log Z(eni)/(V^T*) in [MeV]^^. Results for different with solid lines. Dashed lines: 
Corresponding diagrams from the heavy p limit, see Tab. I and Eq. (25). 

Fig. 16 shows the numerical results. For every diagram, the contribution of the dynamical p-mcson at its physical 
mass of nip ~ 770 MeV (indicated with "770") is displayed. Additionally, the amplitudes for p- masses of = 
1070, 1770, 2770, 10000 MeV are evaluated, multiplying the result with {rn^p/mpf (gray lines). This would correspond 
to a p-meson with mass whose strong coupling g is increased by (rn^p/mp). This is indeed equivalent to the heavy 
p limit from Sec. Ill and convergence of the results from Eq. (C2) towards the heavy p limit of Sec. IV B (dashed 
lines) is observed. This convergence is, on the other hand, a useful tool to check the results from Eq. (C2). 

The large difference of both models at nip — 770 MeV in case of the diagrams (4a) and (eff3) is due to terms that 
partially cancel: diagram (eff3) ~ D{2D — C). For the calculation of the entropy in Eq. (4), the correction to log Z is 
needed which is very different for diagram (eneff) and diagram (enl) as Fig. 16 shows. The discrepancy can be traced 
back to the different high energy behavior of the amplitudes. In any case, the total size of the entropy correction, 
compared to the result of the free pion gas, Eq. (20), is small and of no relevance for the final results. 

2. Calculation of diagram (la2) 

The calculation of one of the diagrams from Fig. 3 is outlined in more detail. The evaluation of the other diagrams 
is carried out in an analog way, with the results given in Eqs. (C2,C4). For diagram (la2) it is most convenient to 
treat the vertex correction first, that is given by the left side of the diagram. The external photon momentum has to 
be set to zero from the beginning of the calculation as has been shown in Appendix A; the matter part of the vertex 
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+ . 



correction reads for an external tt 

oo 1 ioo+e ^ 

rO[fcO - 0, k ^ 0] = fdqq^ f dx — I dq" ^'^ 



—1 — 'ioo+e 



\2 



4m2 - m2 + 2 ((p")^ + (g")2 _ ^2 _ ^/2) _ (^/^2) ((^0)2 _ (^0)2 ^ ^2 _ ^/2)^ 
^ ((pO + ^0)2 _ ^2) ((pO _ qO)2 _ ,^2) (C5) 

where the contour integration method of Ref. [30] is used for tlie summation over Matsubara frequencies. In Eq. 
(C5), cj2 — (jf2 ^ 77^2 g^j-^^ ^2 _ p2 ^ ^2 _ 2pqx + m?p. A problem occurs when closing the integration contour in the 

right q^ half plane: The residue at uj from the double pole of the two pion propagators at the same energy is given by 

Res /(z)|.=. = Im. ^^-^[/(z)(z-..)™] (C6) 

at m = 2. The derivative also applies to the denominator in the second line of Eq. (C5) from the p propagator. The 
integrand exhibits then a divergence of the type 

r[fc° = 0, k ^ 0] - I dq at / = 0. (C7) 

J a- q^ 

The divergence affects only the zero-mode p*^ = 0, but when the pion lines are closed later on. in order to obtain 
diagram (la2), the integrals in Eq. (C5) are not defined any more, and one finds poles of the type l/(a— q^) in the 
three-momentum integration. This infrared divergence, for the external photon at fc = 0, occurs in diagrams that 
contain, besides two or more propagators at the same momentum, an additional propagator as in this case the one of 
the p-meson. 

The complication can be most easily overcome with the introduction of additional parameters according to 

2uj^ d(3^^-' {q^Y - iPu^Y ((^0)2 _ ^2)2' 

and performing the derivative numerically after the three-momentum integration. Still, singularities of the 1/q type 
remain, but they are well-defined by the e-prescription in the g'^-integral of Eq. (C5). We can in this case, as well 
as in all other diagrams from Fig. 3, integrate the angle x = cos(p,q) analytically, thus being left with logarithmic 
singularities, that are easily treated numerically with the help of Eq. (C3). 

It has been checked for all diagrams in Fig. 3 that the poles of the p-meson can be omitted: In Eq. (C5) the 
denominator of the second line from the p-propagator produces two single poles in the right q^ half plane. Taking 
these residues into account in the contour integration leads to deviations of less than 1 % of the result for the vertex 
correction, for all values of (p°,p) and up to temperatures T ^ 200 MeV. Intuitively, this is clear since these poles 
produce a strong Bose-Einstein suppression ^ n[mp] and extra powers of nip in the denominator compared to the 
pion pole. This approximation is made for all results of Eq. (C2). See also Sec. VIA where the approximation is 
again tested. 

The rest of the evaluation of diagram (la2) is straightforward up to the introduction of an additional derivative 
parameter in the same manner as above. As one can see in Fig. 3, a topologically different structure, diagram (Icl), is 
possible for the combination of two ^tttt and two pTTTr-verticcs. This diagram is easily evaluated and has to be added. 



3. The 77rp system at finite ^l. 

Explicit results for logZ from the diagrams (b), (c), and (d) from Fig. 5 are given from which the electric mass 
can be directly calculated using Eq. (18). As argued in the main text, the diagrams (b,c,d) from Fig. 5 lead to the 
same CF as all diagrams with dynamical p from Figs. 3 and 4. For diagram (b), the result is 



X [ — /i] + n_7i[(jj' + ji]) \ogi + {n^n[uj' + /.t] + 7i_n[cj' — /i]) log2 ] (C9) 
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FIG. 17: Elementary diagrams to be saturated with an additional photon in order to construct the selfenergy from Fig. 3. 



with U± and V± from Eq. (27), n± = n[uj ± /i] + 2n[ct'], and 



logi = log 



'ni^p + {p + qY — + uj'Y 

,2 I ^2 



, log2 = log 



ml + {p- qf -{uj- u'f 
'nn?p + {p + qY — — uj'Y 



(CIO) 



with u) ~ q + m^, uj = p + to^, and n the Bose-Einstein distribution. We have checked that logZ(b)(M = 0) = 
logZ((,,ii) from Eq. (C2). The diagram (c) in Fig. 5 with 



iogzmf,) = -^{v+-v.f 



(Cll) 



is zero for fi = and therefore log Z(c) does not contribute to the entropy but only to the CF. 

The diagram (d) in Fig. 5 contains a ppmr vertex that comes from Eq. (9). This interaction is required by the 
gauge invariance of the p with the contribution to md given by 



log Z-[^^{p) = -^^^ [{Ul + U^+ 2DP){Ul + U1 + 2£>^) - 4L"'£>'^] . 
The upper index for U, D indicates which mass has to be used in the meson energy lo of Eqs. (17) and (27). 



(C12) 



4. Charge conservation 



In a calculation of CF the conservation of charge is essential and therefore gauge invariance of the diagrams must 
be ensured. The set of diagrams in Fig. 2 has been constructed using the Ward identity following the procedure 
outlined in Appendix A. They ought to be charge conserving by construction. Nevertheless, it is desirable to have an 
explicit proof. The diagrams from Fig. 2 represent the heavy p limit of the ones with dynamical p-mesons in Fig. 3 
as shown in Appendix C 1. Therefore, it is enough to show charge conservation for the latter. 

From Ref. [63] we utilize the part of the proof that concerns closed loops. The main statement extracted from Ref. 
[63] is, adapted to the current situation: Define a diagram with one external photon at momentum k, not necessarily 
on-shell. By inserting another photon in all possible ways in the diagram, a set of new diagrams of photon selfenergy 
type emerges: For example, the four diagrams in Fig. 17 lead to the photon self energies in the two left columns of 
Fig. 3 plus the (vanishing) diagrams (lc2), (4c), and (5c2) from Fig. 4, once saturated with an additional photon 
(we do not allow direct ^pp and 77pp vertices). The selfenergy diagrams Ilf constructed in this way are charge 
conserving, and fc^; Ilf = for the sum over all diagrams. 

For this statement, it has to be shown first that indeed the diagrams from Fig. 3, including all symmetry and 
isospin factors, turn out from the ones of Fig. 17. This short exercise reveals that there are two classes of selfenergy 
diagrams: one comes from inserting photons in diagrams (1) and (2) of Fig. 17 and the other one from inserting 
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photons in (3) and (4). Thus, there are two separate gauge- invariant classes. In a second step, one has to show the 
statement from Ref. [63] for the current theory which is different from QED and richer in vertices of different type: 

(I) The 77r7r couphngs in Fig. 17 can be transformed into 777r7r couplings by inserting an additional photon. The 
pTTTT vertex can be transformed into a "fpinT vertex. These transformations which are a consequence of the momentum 
dependence of the vertices are essential for the proof. 

(II) For this proof we do not allow direct 7pp and 77pp couplings. However, diagrams which include these couplings 
as in Fig. 4 form a disjoint gauge class anyway. 

(III) The gauge invariance of the diagrams with dynamical p in Fig. 3 survives in the heavy p limit: According to 
Appendix C 1, the amplitudes at a p-mass of are multiplied by (m^/mp)^, with nip the physical mass. Then, the 
limit oo is taken and the effective diagrams of Fig. 2 turn out. The gauge invariance of these diagrams follows. 

This simple graphical proof demonstrates the charge conservation, and, turning the argument around, provides a 
useful tool to ensure that the amplitudes, including symmetry and isospin factors, have been correctly determined in 
Eqs. (C2,C4) and Tab. I. 



APPENDIX D: STRUCTURE OF THE LOW DENSITY EXPANSION 



For a motivation of Eq. (31), we consider the general expression of Eq. (29) for the case of two interacting particles. 
This result is obtained in [42] after carrying out the trace over particle states, i.e. two integrations over the momenta 
Pi and P2 of the interacting particles: 

P f d^pi 1 f d3p2 1 



^^^^ = = - 1 y (# M i (# M ' ' "^^^^'^^^ ^^^^ 

where tjj. 2 = y^Pi 2 + ^'^J- The momenta pi and p2 are defined in the gas rest frame. For simplicity, we consider here 
only a real T-Matrix for the interaction and set /i = 0. The extension to finite fi and complex T is straightforward. 
Note that in the current normalization, T is connected to TJ according to T (327r) TJ (see Eq. (B5)). The 
integrations in Eq. (Dl) can be rewritten in terms of the momentum of the 2-particle cluster, k = pi + p2, and the 
relative momentum in the two-particle cm. frame, Q, which implies a Lorentz boost along k. Using the fact that 



d^pi f d^p2 f d^Q d^k 



2uj[ J 2uj'2 J E^/W+W 



(D2) 



where E = s^^"^ is the total energy of the pious in the cm. system, Eq. (Dl) can be rewritten with a result 
corresponding to Eq. (29). For the moment we ignore the symmetrization operator A in Eq. (29) which will be taken 
care of below. The Lorentz boost of the statistical exponents in Eq. (Dl) is in this case easy to carry out and leads 
to the factor ^^P^'^'^+e^ by noting that the invariant energy is given by E"^ = s = [uj'i + ^2)^ — k^- 

Obviously, no quantum statistical information has entered in Eq. (Dl). However, in Sec. VIIB of Ref. [42] it 
is shown that the fermionic or bosonic nature of the particles can be (partly) included by summing over exchange 
diagrams, i.e., permutating the particles. The final result of this procedure is the replacement of the statistical factors 
e~^" in Eq. (Dl) by Bose-Einstein factors leading to 

u ( n\ ^ f ^'Pi 1 f ^'P2 1 1 1 ^ 

^^(^ = °) = -27 (2^Mi W M^^^^^^^^^'^^^' ^ 

which formally have the appearance of Bose-Einstein factors as shown in Sec. VIIB of [42] (see also the example in 
Sec. VIIC of [42]). As before in the evaluation of Eq. (Dl), we can use at this point Eq. (D2). In order to obtain 
the final form of Eq. (31), finite charge chemical potential and complex T-matrix elements, projected over angular 
momentum, are straightforward introduced. Also, it is more convenient to express the pion scattering in terms of 
isospin amplitudes. The final result is shown in Eq. (31). 

The Lorentz boost into the cm. frame at velocity v = k/(£^^ -I- k^), implicitly contained in Eq. (D2), has also 
to be carried out for the statistical factors in Eq. (D3). Unlike in the simple case of Eq. (Dl), this leads to the 
slightly more complex expressions shown in Eq. (32). Note the boost is not essential but convenient as the scattering 
amplitude is easily obtained in the two-particle cm. frame and some integrations can be carried out analytically. 

Formally, exchange diagrams are included in Eq. (29) through the symmetrization operator A. Note, however, that 
in the standard form of the virial expansion, Eq. (30), effects from exchange diagrams are missing. This highlights 
again the difference between the virial expansion, Eq. (30), and the low density expansion in Eq. (31). 



28 



In fact, Eq. (D3) is not an unfamiliar expression: Let us put T = A =const, i.e. using cj)'^ theory with C = —A/4! 0^, 
and calculate thermodynamic observables such as the pressure from i?2. Using the same interaction, we can compute 
the observables also from thermal loops in the imaginary time formalism (sec e.g. Sec. V) at order A. Results are 
identical. In Sec. VIA this agreement is reconfirmed for more complex interactions. 

The observation of equivalence of the methods from [42] and the thermal loop expansion is, to our best knowledge, 
novel; although in [28], using an effective range expansion for the amplitude, a similar equivalence has been found on 
the level of Eq. (30), i.e., without including Bosc-Einstcin statistics through exchange diagrams. 



APPENDIX E: SOLUTIONS FOR THE RESUMMATIONS 



An additional technical complication appears in the evaluation of Eq. (37) for the summation (n) when the structure 
of the vertices between 7r='=7r='=-loops or 7r°7r^-loops is inspected: The interaction of Eq. (12) leads to a Feynman rule 
of the form (p^ + + 6pq) for the vertex between two charged pion loops of momenta p and q, and of the form 
{p^ + Q^) between a charged and a 7r°-loop, always implying the corresponding shift p° — > p° ± /x (g° q° ^fJ-) for the 
inclusion of finite chemical potential. Therefore, the loops can not be factorized easily in the way Eq. (37) suggests. 
In order to cast the resummations in a manageable from, we introduce for every term of the sum (p^ + q"^ + 6pq) 
an entry in an additional index that runs from 1 to 3. The Eq. (37) is then to be read as a matrix equation in its 
variables. With the definitions 



oo oo 

W± = — / dq UJ n\uj ± /i], X± = — dq ijp'nVo ± //], 



(El) 



additional to the ones of Eq. (17) and (27), the entries of the Faddeev-like equations (37) can be cast in the form 



a± 
Co 

c± 



(Z?, valD, 0) , 

\ {U+ + U-,ml{U+ + U-),V6 (V- - V+)) 




^ ( m2 ([/+ + [/_) 

/ C-3D ml{C-5D) 
J^{C-D) C-3D 
Vo" 



/ W+- 3U+ + W- 



3C/- 



ml {W+ - 5U+ + 5U-) V6{3V+ - X+ - 3V- + X- 



(W+ - C/+ + - C/_) W+ - iU+ + W-- 3C/- 



{X^ -V--X^ 



(E2) 



\^AX--V--x, 



^/Q{:iV+-X+-'iV-+X^) Q{W+ + W^) 



With this extension Eq. (37) is easily solved. In order to check for bulk errors, one can expand the result in the 
coupling constant, and at order /m?p Eq. (26) indeed turns out. At order g'^/rn^p the expansion gives the linear 
chain of three loops which also emerges from the diagram (r) at that order, and the results are identical. 
The ring diagram (r) from Fig. 11 with N "small" loops is given by 



fog^(r),Ar(M) 

-{-1)^I3V 
2Ntt^ 



dp p Res 



n±(p°) 



(pO + fl)^ - 



N 



n[p'] 



n±(-/ 



N 



n[p^] 



no(p") 

(p0)2 - 



N 



n[/] 



(E3) 



for > 3. The residue is taken for the variable at the poles of order N in the right p° half-plane. The tadpole 
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FIG. 18: Resummation (f) from the expansion of the LO chiral Lagrangian to aU orders. 



selfenergies II-i- and Hq for the charged and neutral pion propagator in Eq. (E3) are given by 

n±(p°) = ([(p" + f^f - + 2ml] {U+ + U^+2D]+6 {p" + f,){V. V+)) , 



no(/) = - + 2ml) {U+ + [/_) (E4) 



which is immediately obtained from a± and oq in Eq. (E2). 

There is an additional possible resummation scheme displayed as (f) in Fig. 18. The interaction is obtained from 
the kinetic term of the LO chiral Lagrangian Eq. (6) by expanding it to all orders in the pion fields which means an 
exact calculation of the exponentials U — cxp{i^ / f^) in Eq. (6). The mass correction with A4 from Eq. (6) is tiny 
(see Sec. Ill) and can be safely neglected. The Lagrangian for 2n fields is then given by [n > 2): 

= \2n)\ + 2.+7r-)"-^ ((7r+ d, .-f - 2{.- d, .°)(7r+ 9^' vr")). (E5) 

At /i = 0; the grand canonical partition function from this interaction, summed over all n, results in a surprisingly 
simple expression, 

logZ(f)(^ = 0) = ^ (/^(l - e-t) - D) (E6) 
with D from Eq. (17). This is the special case of the result for finite fi {n > 2), 



4D(n - fc - 2)! iml{U+ + U-f - (2 + fc - n){V+ - V-f) 
T{-\/2-k){U- + U+) 

r(n - k){^ml{U+ + U^){2D + U+ + [/_) + (1 + fc - n){V+ - V^f) -T{l-k + n){V+ - V^)^ 

r(i/2-fc) 



(E7) 



The sum over fc comes from the expansion of the polynomial of order n — 2 in Eq. (E5). The possibilities of contracting 
2fc neutral pion fields have been rewritten, Hi^o k-iC^^ - 1 - 2i) = (-2)'=r(l/2)/r(l/2 - fc). 

The structure of the Lagrangian in Eq. (E5) resembles the vertex structure of resummation (t) from Sec. VII 
with the ppinr interaction: At any order n > 2 in the interaction, there are only two derivative couplings. Also, the 
diagrammatic representation of resummation (t) has the same topology as diagram (f) once the heavy p limit is taken. 
Indeed, we observe a close numerical correspondence between the resummations (t) and (f). Thus, it is interesting to 

(2) 

note that the p tadpole resummation is well described by an expansion of £W to all orders. The resummation (f) is 
not included in the final numerical results due to these potential double counting problems with (t). 
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APPENDIX F: EXTENSION TO SU{S) 

It is straightforward to extend the study of CF and other thermodynamical observables to SU{3). Compared to the 
pion, the other members of the meson octet have higher masses which simphfies the selection of relevant processes in a 
thermal heat bath: we regard diagrams which do not contain any pion as kinematically suppressed. The contribution 
to log Z at then consists of diagram (b) in Fig. 5 with one pion line replaced by a kaon and the p replaced by the 
K*{892). The ttKK* interaction follows from Eq. (8) in the SU{3) version by extending the representation in Eq. 
(7) to the full meson and vector meson octet in the standard way [34, 35]. The result reads 



logZ7,f(M) 



9^f3V 
32 



128^4 l^™- 



2m 



K 







dp / dq 



pq 



+ ( n+n[(jj' + /.t] + n_n[ti)' 



fj] + n^n[uj' + m] + 2 ^[^'] {''^[^ + m] + '^['^ ^ m]) ) logi 



(Fl) 



where lo"^ = q^ 



m 



,'2 



P 



"^^j ^± = Ti[(X> ± n] + n[uj], and the upper index specifies the mass that has to be 



used in the definitions of D and U from Eqs. (17) and (27). The expressions logj^ and log2 are given by Eq. (CIO) 
with the replacement nip itik- and u, u' defined as in Eq. (Fl). 

Diagram (c) from Fig. 5 with tt, p, and K is possible. Also, the K* K*t:t: term from Eq. (9) is present, shown in 
Fig. 5 (d) with the p replaced by a i4'*(892). The corresponding contributions are 



8m 

3gW r 
32 



+ U^' + 2D'^' 



(F2) 



The electric mass from Eqs. (Fl) and (F2) is plotted as "ttK dynamical" in Fig. 13. 

The ttK interaction can be alternatively described by the LO chiral Lagrangian from Eq. (6) in the SU (3) version 
(we do not try to construct an effective, point-like, ttK interaction from K*{892) exchange as it has been done for 
TTTT via p exchange). Using similar arguments as above, the calculation is reduced to diagram (a) in Fig. 5, with one 
pion replaced by a kaon. Taking only the kinetic part of Eq. (6) — contributions from the mass term are tiny — one 
obtains 



~f3V 



logZ(-f (m) = —f-r- [6 {V; V^) {Vl' Vl') + {ml + ml) {U^' + [/^' + 2D^) {U^ + Ul + D^)] 

"O/tt JK 



(F3) 



with /if = 1.22/^ taken from chiral perturbation theory [64]. The contribution from Eq. (F3) is plotted as ^'ttK 
contact" in Fig. 13 with the dotted line. 

In a similar way as in Sec. VI, it is possible to establish a density expansion for the ttK interaction that respects the 
Bose-Einstein statistics of the asymptotic states in ttK scattering. Following the same steps as in Sec. VI, we obtain, 

1/2 

again assuming elastic unitarity (the 50 % inelasticity in the 62 partial wave changes the result only slightly). 



B. 



{ttK), Bo 



= 4^ 



dE / dx / dk 



E (2£+l) 



'=0,1,2, 



X 2 



S^^ {n[ujTr + ^J]'>^■[uJK + m] + n-i'^TT - p]n[uJK -~ p]) 
+ ^ (sy^ + 2(5^/^^ {n[uj^]n[ujK + A«] + n[Lj^]n[Lj/f - p]) 
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/2 _L a3/2 



{n[uj^ - p]n[ujK + m] + n[^T: + lAn[uJK - p]) 



(F4) 
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The boosted Bose-Einstein factors are 

1 / kQx \ ( kQx 



7/ = ( 1 - ;g2^ ) ' E-^K = + = 2^ (F5) 



with the cm. momentum of the particles, Q ~ \/{2E)^{E'^ — (m^r + mK)^){E^ — [m-jr — niK)^)- For fj, = and in 
the Bohzmann hmit Eq. (F4) reduces to the virial coefficient 



Tj{-^K), Boltz 



OO 

(/^ = 0) = ^ / dEE^K.iPE)^ (2^+l)(4<5f' + 25y') (F6) 

^ /I — n 1 n 



m^+m/f £-0,1,2 



3/2 1 /2 

which shows the correct ratio of degeneracy between S/ and 6^ , but is an overall factor of 4 larger than one would 
expect — compare, e.g., to Eq. (30): Instead of a sum over isospin of the form J^i ^(2^+ 1)(2^+ l)6j^, the projection 
of charge channels of pions and kaons to the isospin channels leads to ^(2/+ 1){2£+ 1)6^. In any case, the result 
Eq. (F4) for TOcI, using the chiral ttK interaction at l/if-TrfK), matches exactly the thermal loops in Eq. (F3). This we 
have shown in the same way as in Sec. VIA by using Eq. (35) and the partial waves T^/^ = {7u — 5s — 2t) / [12 f f k) 
and T^/'^ = (2s — t — u)/{6fT^fK) (as in Eq. (F3), we consider only the kinetic term of d'^x)- A similar test has 
been performed by starting from Eq. (F6) and using the partial waves from above. From this the pressure has been 
calculated and results are identical to the pressure obtained from Eq. (F3) by taking the Boltzmann limit of the 
statistical factors n in the definition of D, J7, and V . Additionally, an independent check for the Lorcntz structure of 
Eq. (F4) has been performed in the same way as in Sec. VIA, this time for a (/)i02 interaction of uncharged bosons 
with different masses m^^ and m^^. 
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